/* 
 * File:   DinamicaMolecular.cpp
 * Author: jonk
 *
 * Created on 15 de enero de 2012, 20:19
 */

#include <cstdlib>
#include <fstream>
#include <iostream>
#include <string.h>
#include <iomanip>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include<time.h>
#include<iomanip>
#include <sstream>

#define pi 3.1415926535898
#define np 200
//#define CONFINED
using namespace std;

#include "objetos.h"
#include "mthParticulas.h"
#include "ficheros.h"
#include "inicial.h"
#include "inter.h"
#include "integracion.h"

int main(int argc, char** argv) 
{
    parametros param;
    particula *part = new particula[np];
    double epoten;
    param.u0=100.0;
    param.h0=1.0;
    param.omega=0.01;
    param.phi=0.1;
    param.rcut=2.0;
    param.t=0.0;
    param.nacc=0.0;
    param.zeta=10000.0;
    param.dt=0.00005;
    param.npart=np;
    
    /*Colocamos las particulas*/
    SetRandomPart(part, &param);
    
    for(int i=0;i<100000;i++)
    {
        param.nacc=i;
        param.t=i*param.dt;
        epoten = integracion(part,  param);
        if(i%10==0)
        {
                cout<<epoten<<endl;
                printConfigBackup(part, param);
        }
        
    }
    
    
    delete []part;
    return 0;
}

